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Abstract 

In the present paper, we introduce a numerical scheme for the price 
of a barrier option when the price of the underlying follows a diffusion 
process. The numerical scheme is based on an extension of a static 
hedging formula of barrier options. For getting the static hedging for- 
mula, the underlying process needs to have a symmetry. We introduce 
a way to "symmetrize" a given diffusion process. Then the pricing of 
a barrier option is reduced to that of plain options under the sym- 
metrized process. To show how our symmetrization scheme works, 
we will present some numerical results applying (path-independent) 
Euler-Maruyama approximation to our scheme, comparing them with 
the path-dependent Euler-Maruyama scheme when the model is of the 
Black-Scholes, CEV, Heston, and (A)-SABR, respectively. The results 
show the effectiveness of our scheme. 



Keywords: Numerical Scheme for the Barrier Options, Put-Call 
Symmetry, Static Hedging, Stochastic Volatility Models. 

1 Introduction 

In financial practice, the pricing (and hedging) of barrier type deriva- 
tives becomes more and more important. In the Black-Scholes envi- 
ronment, some analytic formulas are available in iMerton ()l973h . 
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If the underlying process is a diffusion process which is more com- 
plicated than a Geometric Brownian Motion, it will not be able, basi- 
cally, to rely on anymore an analytic formula. Instead one should 
resort to some numerical analysis. There is a problem, h owever, 



arisin g from the path-dependence of the pay-off function. iGobet 



(2000) pointed out that the weak convergence rate against the time- 



discretization gets worse compared with the standard path-independent 
pay-off cases due to the failure in the observation of hitting between 
two time steps. He showed that the weak order of Euler-Maruyama 
approximation is \ , which is much slower than the standard case where 
the order is 1. 

In the present paper, we will introduce a new numerical scheme 
where the pricing (and hedging) of barrier options are reduced to that 
of plain (pat h- independent ) ones . The scheme is based on an observa- 
tion made by Carr and Lee! ( 20091 ) which we will refer to as "arithmetic 



put-call symmetry" (APCS). In the Black-Scholes economy, it is well- 
reco gnized that the r eduction is possible due to the reflection principle 
(see Imamura ( 201ll )). The put-call symmetry is an extension of the 



reflection principle, with which a semi-static hedge is still possible. 
There are two keys in our scheme; 

1. For a given diffusion X and a real number K, we can find another 
diffusion X which satisfies the APCS at K (see section |23|) . We 
call this procedure "symmetrization" . 

2. For T > 0, the expectations E[f(X T )l {T>T} ] and E[f(X T )l {f>T} ] 
coincides, where r and f are the first hitting time at K of X and 
X, respectively. 

We do not anymore regard the equation for semi-static hedging but 
just a relation to calculate the expectation of the diffusion with a 
Dirichlet boundary condition in terms of those without boundary con- 
ditions. In other words, the pricing and hedging are reduced to path- 
independent ones, where many stable techniques are available. In 
this paper, we will present some numerical results of applying (path- 
independent) Euler-Maruyama (EM) approximation to our scheme, 
comparing them with the path-de pend e nt EM under Constant Elas- 
ticity of Volatility (CEV) models (|Coxl ^197$ )) including as a special 



case the Bla ck-Schq l es (B S) model, and stoc hastic volatility models 



(|2002! ): lHenrv-Laborderel (120051 )). 



of He st on's (IHestonl (11993 ) ) an d (A)-SABR (jHagan and Woodward 



This paper consists of two parts. In the first part, the discussion of 
our new scheme is concentrated on one-dimensional diffusion models, 
while the latter part deals with applications to the stochastic volatil- 
ity models. Mathematically, the first part is somehow self-contained, 
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w hile one may think the latter part to be dependent on the result 



m lAkahori an d Imamura (2012). The fa ct is that we have found, in 
advance of Akahori and Imamura ( 20121 ). through numerical experi 



ments how it should be applied to stochas t ic vo latility models (see 
Imamura. Ishieaki. Kawaeoe. and Okumural (|2012l )). In anyway, the 



main aim of the present paper is to introduce the new scheme and to 
report numerical results which show the effectiveness of the scheme. In 
order to ensure the consistency of the experiments, we present detailed 
descriptions. 

The paper is organized as follows. In Section [2l we recall the 
APCS and how it is applied to the pricing and semi-static hedging 
of barrier option s . In s ection 12.21 we give a sufficient condition shown 
bv lCarr and~Lel (|2009l ) under which APCS holds. In Section [2.31 we 



introduce a way to "symmetrize" a given diffusion process. We then 
show that by using the symmetrized process which satisfies APCS, 
the pricing of barrier option is reduced to that of two plain options. 
In section we give numerical examples under our symmetrized ap- 
proximation method. The results of the path-wise EM scheme (in 
section IBTTj) and our new scheme are given when the underlying asset 
price process follows CEV with the volatility elasticity f3 = 1 which 
is nothing but a BS model and other elasticities. From Section U we 
discuss applications to stochastic volatility models. The symmetrized 
method is also applicable to the stochastic volatility models where 
the underlying price process and its volatility follows a (degenerate) 
2-dimensional diffusion process. In section [4~Tl we give numerical re- 
sults under Heston model and A-SABR model. In section 14.21 we 
show that the symmetrization scheme also works for the pricing of 
double barrier option. One will find that our scheme overwhelms the 
path-wise EM in all numerical results. 
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2 The Put-Call Symmetry Method for 
One Dimensional Diffusions 



2.1 Arithmetic Put-Call Symmetry 

Let X be a real valued diffusion process defined on a complete filtered 
probability space (f2, T, P, {.Ft}) which satisfies the usual conditions. 
For fixed K > 0, we say that arithmetic put-call symmetry (APCS) at 
K holds for X if the following equation is satisfied ; 

E[G(X t - K) | X = K] = E[G(K - X t ) \ X = K], 

for any bounded measurable function G and t > 0. The APCS at K is 
alternatively defined to be the equivalence in law between X t — X t/KT 
and Xt/\ T — Xt for any t > 0, where r is the first hitting time of X to 
K. 

Intuitively, the APCS means the following. For every path of X 
which crosses the level K and is found at time t at a point below K, 
there is a "shadow path" obtained from the reflection with respect to 
the level K which exceeds this level at time t, and these two paths have 
the same probability. For one-dimensional Brownian motion, APCS 
holds for any K > since the reflection principle holds. 

In ICarr and Lei l[200flh . the APCS, or more precisely, PCffj, is 



applied to the pricing and semi-static hedging of barrier options. Semi- 
static hedging means replication of the barrier contract by trading 
European-style claims at no more than two times after inception. 

In more detail, we have the following; if X satisfies APCS at K, 
then for any bounded measurable / and T > 0, 

E[f{X T )I {r>T} } = E[f(X T )I{x T>K> r>T }] 

= E[f(X T )I{x T >K}] - E[f(X T )I{x T >K, t<T}\, 

where 

r = inf{t >0:X T <K}. (1) 
By APCS of X, we see that 

^[f( X T)I{X T >K, r<T}} = E [ E [/(^r)^{X T >A'}l- ;r T]^{T<T}] 

= E[E[/(2K - X t )I {Xt<k] \F t \I {t < t} \. 
Hence we obtain the following equation; 

E[f(X T )I {T>T} ] = E[f(X T )I {XT>K} ] 

-E[f(2K-X T )I {XT<K} }. 



1 They defined PCS as the equality of the distribution between ^ under P and ^ 
under Q, where ^ = 
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Of the equation ([2]), the left-hand-side reads the price of a barrier 
option written on X, whose pay-off is /, knocked out at K, and the 
right-hand-side is the price of a combination of two plain- vanilla op- 
tions. 

Here is a description of the hedging strategy of a barrier option 
implied from the right-hand-side of 

1. Hold a plain- vanilla options whose pay-off is /(Xt) if the price 
at the maturity is less than K, and is nothing if the price at the 
maturity is greater than K. Moreover short-sell a plain- vanilla 
options whose pay-off is f(2K — Xt) if the price at the maturity 
is greater than K, and is nothing if the price at the maturity is 
less than K. 

2. Keep the above position until the price hits the barrier K. If the 
price never hits K until the maturity, the pay-off is /(Xt)- 

3. If the price hits K, clear both plain-vanilla positions at the hit- 
ting time. Indeed, the value of two options are exactly the same 
at r. 



2.2 APCS of diffusion process 

Let X be a solution to the following one-dimensional stochastic dif- 
ferential equation (SDE) driven by a Brownian motion W in R, 

dX t = a(X t )dW t + fi(X t )dt. (3) 

Here we assume the following hypotheses; 

(HI) a : R — > R and fj, : R — > R are locally bounded measurable 
functions such that the linear growth condition is satisfied, ie, 
for a constant C, \cr(x)\ + \^(x)\ < C(l + |x|) for any x £ R. 

(H2) The following condition is satisfied; 

a (y) 7^ < ^ = ^ > <T ~ 2 is integrable in a neighborhood of y. 

Then we have the following result on the uniqueness of the solution 
to ©; 

Theorem 2.1 (Theorem 4, Engelbert and Schmidt ( 19851 )). Under 



(H2), there exists a unique (in law) solution satisfying SDE 

Moreover, by the linear growth condition (HI), the unique (in 
la w) solution will not ex plode in finite time. 

Carr and Lee (2009) gave a sufficient condition for a solution of 



([3D to satisfy PCS at € R. The following Propo sition is essentially 
a corollary to Theorem 3.1 in ICarr and Leel l[2009l ). 
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Proposition 2.2. If the coefficients further satisfy the following con- 
ditions; 

a(x) = e{x)a{2K - x) (x e R\{if}), (4) 
for a measurable e : R — > {— 1, 1} and 

H(x) = -n(2K-x) (x G R\{if}), (5) 

i/ien ylPCS at K holds for X. 

Proof. By the uniqueness in law, it is sufficient to show that {Xt/\r ~ 
{X t — X t /\ T ))t>o solves the SDE ([3]). By the assumptions dH) and (JSJ), 
we obtain that 

X tAT - (X t - X tAT ) = X tAT - f a(X s )dW s - f fi(X s )ds 

JtAr JtAr 

= X tAT - f e(X s )a(2K - X s )dW s 

JtAr 

+ / n(2K-X s )ds 

JtAr 

= X tAr - f e{X s )a(X T - (X s - X T ))dW s 

JtAr 

+ ! n(X T - {X 9 - X T ))ds. 

JtAr 

We set W( = W t Ar — !tAT £ (Xa)dW s . Since we obtain that 

(W')(t) = (W)(t) 
= t, 



W' is a Brownian motion (cf. Ikeda and Watanabe ( 19891 ) Chapter II 



Theorem 6.1.). Therefore we see that 



XtAr - {X t - X t Ar) = X + a{X sAr - (X s - X sAT ))dW' s 

Jo 

+ / n(X sAr - (X s - X sAr ))ds. 
Jo 

Hence APCS at K holds. □ 

2.3 Symmetrization of Diffusion Processes 

We introduce a way to "symmetrize" a given diffusion to satisfy APCS. 
By using this symmetrized process satisfying APCS, the pricing of 
barrier options is reduced to that of plain options. 
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We start with a diffusion process X given as a unique solution of 
SDE ([3]). We do not assume that the coefficients have the symmetric 
conditions (JU) and ©. We then construct another diffusion X that 
satisfies APCS at K in the following way. Put 

I a{2K — x) x < K, 

~ (x ) . = \^ x ) x > K 

^ [X> ' \-n(2K-x) x < K, U 

and consider the following SDE; 

dX t = a(X t )dW t + fi(X t )dt. (8) 

Again by Theorem 12. 1[ there is a unique (in law) solution Xf. Then 
we obtain the following result; 

Theorem 2.3. It holds that 

B[f(X T )I {T>T} ] = E[f(X T )I { x T>K} ] 

-V{f(2K-X T )I {jtT<K} ). 

Proof. Since a and jl satisfy the condition ^ and APCS at K 
holds for X by Proposition 12.21 Then the equation ([2]) is valid for 
X. Moreover, by the definition of a and jl, we have a(x) = a(x) and 
fj,(x) = fi(x) for x < K. Therefore by the uniqueness in law of the 
SDE, we have that {X t }t< T = {Xt}t<r pathwisely. Then we see that 
t = f where f = inf{t > : X% < K}. Hence we have 

E[/(Xr)/ {f>T} ] = E[/(X T )7 {T>T} ]. 

□ 



2.4 Important Remark 

We do not anymore regard ([9]) as an equation for semi-static hedging 
but a relation to give a numerical scheme to calculate the expectation 
of the diffusion with a Dirichlet boundary condition in terms of those 
without boundary conditions. The former is very difficult while the 
latter is rather easier using rapidly developing technique from numer- 
ical finance. In the following sections, we will present some results of 
numerical examples to show the effectivity of the new scheme. 
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3 Numerical Experiments for One Di- 
mensional Models 



3.1 The Euler-Maruyama Scheme 

Here we briefly recall the Euler-Maruyama scheme for a diffusion pro- 
cess given as a solution to SDE ([3]). Fix T > 0. For n > 1, we set a 
subdivision of the interval [0, T] 

= t < h < ■ ■ ■ < t n = T, 

where tj. := ^ for < k < n, and we denote this net by A n . 

The Euler-Maruyama scheme is a general method for numerically 
solving © by a discretized stochastic process which is given by 

K+i = x l + - **) + - Wt„), (io) 

k = 0, 1, 2, • • • , n — 1, and for i& < i < t^+ii -^t is given by an 
interpolation. The approximating process (X^) is simulated by us- 
ing independent quasi-random Gaussian variables for the increments 

(Wtfc+i ~ w t k )o<k<n-i- 

We rely on the following result; 



Theorem 3.1 (Theorem 3.1. lYanl ((2002)). If the set of discontinuous 



points of a and fi is countable, then the Euler scheme ITOj) converges 
weakly to the unique weak solution of SDE (0) as n — > oo. 

From now on, in addition to (HI) and (H2), we assume that a 
and fi are piece-wise continuous. 



3.1.1 Path- wise Method 

Since the convergence is in the space of probability measures on con- 
tinuous functions, we see that this algorithm can also be used to 
simulate a path-dependent functional of the process; in particular, 
f{XT)I{ T >T}-> where / is a (bounded) continuous function and r is 
the first hitting time defined by ([1]). The functional is approximated 
by f(Xj,)Ir T n >T \, where r n := inf{£fc : < K} is the discretized 
first hitting time to K. Then the expectation Hj[f(XT)I{ T >T}] is a P" 
proximated with a Monte-Carlo algorithm by 

Method 1. (Path-wise EM scheme) 
1 M 

-J2f( X T^i))hr^i)>T}- (11) 
i=l 
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By the strong law of large numbers, (fTTj) converges to E[f(X^)I^ T n >T y] 
as M goes to infinity. Moreover, as the index n of the net A n goes to 
in finity, E\££X± )I{T n >T}] converges to E[/(Xt)J{ t >t}]- s According 



to 



Gobetl (|2000l ). the following convergence rate was given; 



Theorem 3.2 (Theorem 2.3. lGobetl (|2000l )). Assume that a and [i are 



in C£°, a is bounded below from zero and a solution is non-explosion. 
Then for a bounded measurable function f such that d(supp/, K) > 0, 
there is a constant C such that 

|E[/(X«)/ {r „ >T} ] -E[/(Xr)/ {r>T} ]| < C-J=- 
3.1.2 Put-Call Symmetry Method 

Let I be a solution with coefficients a and fi given by ([6]) and (J7J), 
and (X™) be the discretized Euler-Maruyama process with respect to 
the net A n . Namely, 



XI 



XI + {a{Xl){t k+l - t k ) + n{Xl){W tk+1 - W tk ))I {jt n >K} 
+ (a(2K - X? k )(t - t k ) - »(2K - X? h )(W t - W tk ))I { ^ K} 

for k = 0, 1, 2, • • • , n — 1. With an interpolation, X™ for t k < t < ifc+i 
is obtained as well. Since the set of the discontinuous in the coefficients 
is Lebesgue measure zero, X n also converges weakly to X by Theorem 

ED 

Combining Theorem l2.3l and l3-H we may rely on the following algo- 
rithm; the expectation E[/(Xt)/{ t>T |] is approximated with Monte- 
Carlo algorithm by 

Method 2. (Put-Call symmetry method) 
1 M 

M 2 {f^T{Ui))I { xn {uJl ) >K } ~ f( 2K ~ ^t(^)K{X ? (^)<E-}} • 

1=1 

(12) 

As M goes to oo, (fT2|) converges to 

E[f(X2)I {jt n >K} ] - E[/(2tf - X2)I { ^< K} }. (13) 
By the weak convergence of X n , (|13p converges to 



E[f(X T )I {jcT>K} ] - B[f(2K - X t )I {Xt ^ k} ] 
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as n — > oo. However, we don't know the exact rate of convergence in 
this algorithm since the coefficients are inevitably non-smooth at 

The numerical results in the next section, however, may imply that 
the convergence rate of Put-Call symmetry method is better than that 
of the path-wise EM scheme. To prove this conjecture would be a very 
interesting mathematical challenge. 



3.2 Numerical Results 

In this section, we give numerical examples using method 1 (path- 
wise EM method) and method 2 (Put-Call symmetry method) under 
Black-Scholes model and other CEV models. Let us consider the value 
of barrier call option with strike price S and knockout barrier K. 



3.2.1 Black-Scholes Model 

The underlying price process of Black-Scholes model is the unique 
solution of the following SDE; 

dX t = rX t dt + aX t dW t , (14) 

for r, a > 0. Then the value of barrier option is accurately-calculable 
since the joint distribution of Brownian motion and the hitting time 
of Brownian motion to a point is computable by using the reflection 
principle. The exact option price is given by the following; 

E[(X T - S)+I {T>T} ] = V call (X ) - ^ V call (^) , 

where 

Vcauix) = xe rT (1 - *(£*+(*))) - S (1 - $(£*_(*))) , 

and $ is the distribution function of the standard normal distribution. 

Fix a maturity time T > 0. Tables [T]d] give simulation results for 
the value of down-and-out call option E[(Xr — 5') + ^{r>T}] under the 
path-wise Euler-Maruyama method (EM) and the Put-Call symmetry 
method (PCM). We take [X Q = 100, S = 95, K = 90, T = 1], and 



2 There are many results on the rate of convergence when a and /i are smooth. For 
example, when a and \x are in Cp, the space of functions such that 4-th derivative exists 
and have a polynomial growth, we have 

|E[ 5 (X£)]-E[ 5 (X r )]|=0(n), 

for any jed (See Kloeden and Platen iKloeden and Plater] ( 2011 ). pp. 476). 
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Table Q] a = 0.2, r = 0, 

Table [2 a = 0.2, r = 0.02, 

Table [3J cr = 0.5, r = 0, 

Table H a = 0.5, r = 0.02. 

In the PCM, we symmetrize the functions fJ-(x) = rx and cr(x) = ax 
at K. 

The number of simulation trials is set equal to the cube of the 
number of time steps for the Euler discretization. The errors in the 
last two columns are calculated as 

|EM(PCM) - true option price] 
true option price 

One sees that, in the experiments, the Put-Call symmetry method 
always beats the path-wise EM method. 

3.2.2 CEV Model 

Here the underlying price process is a solution of the following SDE; 

dX t = rX t dt + aX^dWt, (15) 

for r, a > and /3 > \. We take 0.75 in the experiments. 

Tables 0E1 ar e simulation results for down-and-out call options 
with EM and PCM. We set parameters to [X = 100, S = 95, K = 
90, p = 0.75, a = 0.45, T = 1], and 

Table [5j r = 0, 

Table E r = 0.02, 

in the experiments. For CEV model, we do not have an analytic 
formula. So, as a benchmark, we used numerical results by the path- 
wise Euler-Maruyama scheme where the number of time steps for the 
Euler discretization is 5, 000 and that of a Monte-Carlo simulation is 
50, 000, 000. Note that since we are calculating the prices of down- 
and-out call options, we do not need to care about the singularity at 
x = in the SDE. 



4 Put-Call Symmetry Method Applied 
to Stochastic Volatility Models 

In this section, we slightly extend the put-call symmetry method to 
apply it to stochastic volatility models which are described by two- 
di mensional SDE. Theor e tical backgrounds of the extension is given 
in Akahori and Imamura ( 20121 ). 
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Table 1: Black-Scholes model 

X Q = 100, S = 95, K = 90, cr = 0.2, r = 0, T = 1, option price = 8.17140 



No. of 


No. of 






EM 


PCM 


simulation trials 


time steps 


EM 


PCM 


error (%) 


error(%) 


1000 


10 


8.881 


7.816 


8.7 


4.4 


8000 


20 


9.183 


8.172 


12.4 


0.0 


27000 


30 


8.992 


8.250 


10.0 


1.0 


64000 


40 


8.880 


8.175 


8.7 


0.0 


125000 


50 


8.804 


8.190 


7.7 


0.2 


216000 


60 


8.692 


8.137 


6.4 


0.4 


343000 


70 


8.697 


8.127 


6.4 


0.5 


512000 


80 


8.671 


8.171 


6.1 


0.0 


729000 


90 


8.672 


8.207 


6.1 


0.4 


1000000 


100 


8.597 


8.135 


5.2 


0.4 



Table 2: Black-Scholes model 

X Q = 100, S = 95, K = 90, cr = 0.2, r = 0.02, T = 1, option price = 9.31138 



No. of 


No. of 






EM 


PCM 


simulation trials 


time steps 


EM 


PCM 


error (%) 


error(%) 


1000 


10 


10.953 


9.821 


17.6 


5.5 


8000 


20 


10.050 


9.165 


7.9 


1.6 


27000 


30 


10.090 


9.226 


8.4 


0.9 


64000 


40 


9.952 


9.258 


6.9 


0.6 


125000 


50 


9.974 


9.302 


7.1 


0.1 


216000 


60 


10.033 


9.389 


7.7 


0.8 


343000 


70 


9.911 


9.298 


6.4 


0.1 


512000 


80 


9.885 


9.353 


6.2 


0.4 


729000 


90 


9.839 


9.306 


5.7 


0.1 


1000000 


100 


9.811 


9.309 


5.4 


0.0 
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Table 3: Black-Scholes model 

X Q = 100, S = 95, K = 90, cr = 0.5, r = 0, T = 1, option price = 9.37170 



No. of 


No. of 






EM 


PCM 


simulation trials 


time steps 


EM 


PCM 


error (%) 


error (%) 


1000 


10 


15.981 


9.521 


70.5 


1.6 


8000 


20 


14.455 


9.742 


54.2 


4.0 


27000 


30 


13.074 


9.126 


39.5 


2.6 


64000 


40 


12.837 


9.479 


37.0 


1.1 


125000 


50 


12.281 


9.251 


31.0 


1.3 


216000 


60 


11.942 


9.231 


27.4 


1.5 


343000 


70 


11.838 


9.307 


26.3 


0.7 


512000 


80 


11.750 


9.450 


25.4 


0.8 


729000 


90 


11.549 


9.392 


23.2 


0.2 


1000000 


100 


11.443 


9.319 


22.1 


0.6 



Table 4: Black-Scholes model 



= 100, S = 95, K = 


90, a = 0.5, r 


= 0.02, T 


= 1, option price = 10 


.02470 


No. of 


No. of 






EM 


PCM 


simulation trials 


time steps 


EM 


PCM 


error (%) 


error (%) 


1000 


10 


15.488 


9.688 


54.5 


3.4 


8000 


20 


14.687 


9.540 


46.5 


4.8 


27000 


30 


14.065 


10.341 


40.3 


3.2 


64000 


40 


13.472 


10.191 


34.4 


1.7 


125000 


50 


13.012 


9.779 


29.8 


2.4 


216000 


60 


12.981 


10.257 


29.5 


2.3 


343000 


70 


12.707 


9.991 


26.8 


0.3 


512000 


80 


12.391 


9.916 


23.6 


1.1 


729000 


90 


12.418 


10.098 


23.9 


0.7 


1000000 


100 


12.235 


10.025 


22.0 


0.0 
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Table 5: CEV model 

X Q = 100, S = 95, K = 90, /3 = 0.75, a = 0.45, r = 0, T = 1, 
benchmark of option price = 7.50095 



No. of 


No. of 






EM 


PCM 


simulation trials 


time steps 


EM 


PCM 


error (%) 


error(%) 


1000 


10 


7.781 


7.068 


3.8 


5.7 


8000 


20 


7.997 


7.504 


6.6 


0.1 


27000 


30 


7.805 


7.397 


4.1 


1.4 


64000 


40 


7.758 


7.379 


3.5 


1.6 


125000 


50 


7.730 


7.412 


3.1 


1.2 


216000 


60 


7.733 


7.407 


3.1 


1.2 


343000 


70 


7.714 


7.422 


2.9 


1.0 


512000 


80 


7.691 


7.423 


2.6 


1.0 


729000 


90 


7.680 


7.414 


2.4 


1.1 


1000000 


100 


7.654 


7.414 


2.1 


1.1 



Table 6: CEV model 

X = 100, S = 95, K = 90, /3 = 0.75, cr = 0.45, r = 0.02, T = 1, 
benchmark of option price = 8.82718 



No. of 


No. of 






EM 


PCM 


simulation trials 


time steps 


EM 


PCM 


error (%) 


error (%) 


1000 


10 


9.418 


8.918 


6.7 


1.0 


8000 


20 


9.349 


8.986 


5.9 


1.8 


27000 


30 


9.242 


8.791 


4.7 


0.4 


64000 


40 


9.193 


8.772 


4.1 


0.6 


125000 


50 


9.109 


8.760 


3.2 


0.8 


216000 


60 


9.089 


8.751 


3.0 


0.9 


343000 


70 


9.063 


8.742 


2.7 


1.0 


512000 


80 


9.009 


8.722 


2.1 


1.2 


729000 


90 


9.027 


8.745 


2.3 


0.9 


1000000 


100 


8.995 


8.722 


1.9 


1.2 
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A generic stochastic volatility model is given as follows; 

dX t = au(X t , V t )dW t + M i(X t , V t ) dt 
dV t = a 21 (V t )dW t + a 22 {Vt)dB t + /x 2 (F t ) dt, ' 01 

where W and B are mutually independent (1-dim) Wiener processes, 

a(x v) = (°n(x,v) 

and /J,(x,v) = (fii(x,v), ^{v)) are continuous functions on R 2 . Here 
we simply assume that a and \x are sufficiently regular (not so irreg- 
ular) to allow a unique weak solution in ()16p . The independence of 
V against S plays an important role in applying our scheme. In fact, 
thanks to the property, we may simply work on the symmetrization 
with respect to the reflection (x,y) \-t (2K — x,y). Let us be more 
precise. Let (X, V) be a 2-dimensional diffusion process given as a 
(weak) unique solution of SDE (fT6|) . and r be the first passage time 
of X to K. We note that r is not dependent on V . We say that 
arithmetic put-call symmetry at K holds for (X, V) if 



(X t , V t )l {T < t} ± (2K - X t , V t )l {T < t} 



for any t > 0. 



Mathematically, we rely on the following result from lAkahori and Imamura 



(2012). 



Proposition 4.1 ( Akahori and Imamura ( 20121 )). If the coefficients 



satisfy the following conditions; 

a n (x,v) = -a u (2K - x,v), (17) 

Hi{x,v) = -fi(2K-x,v), (18) 
for (x, v) e (R \ {K}) x R, then APCS at K holds for (X, V). 

On the basis of Proposition 14. 1[ we construct another diffusion 
(X ,V) that satisfies APCS at K in a totally similar way as the one 
dimensional case, and we obtain a static hedging formula correspond- 
ing to Theorem 12.31 

Proposition 4.2. Let K > and put 



a n (x,v) 



Hi{x,v) 



CTii(a;, v) x > K 

-o-u(2K - x,v) x < K 

/xi(x, v ) x > K 

-Hi(2K - x,v) x < K ' 
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and let X be the unique (weak) solution to 

dX t = a n {X t , Vt)dW t + fitiXt, V t ) dt, 

where V is the solution to h!6}) . Then, it holds for any bounded Borel 
function f and t > that 

E [f( X t)l{Xt>K}l{T K >t}\ 

= E[f(X t )l {jCt>K} ] - E[f(2K - X t )l {jit<K} ], 1 •' 



where X is the solution to [To]) with Xq > K . 

Proof. Omitted. □ 

4.1 Numerical Results on Single Barrier Op- 
tions under Stochastic Volatility Models 

In this section we give numerical examples of the price of a single 
barrier option under Heston's and SABR type stochastic volatility 
models, using numerical method based on ([25]) . 

The Euler-Maruyama scheme of the solution of SDE (|16p with 
respect to the net A n = {to,t\, • • • , t n } is given by the following; 

K +1 = x l + <^(K> K)( w t k+ i - w t k ) + mPOL, v£)(t*+i - tk), 

V t n k+1 = Vr k + a 21 (Vr k )(W tk+1 -W tk ) + a 22 (Vr k )(B tk+1 - B tk ) 
+ fM 2 (V t n k )(t k+l -t k ), 

for k = 0, 1, 2, • • ■ ,n—l. With an interpolation, Xt for t}~ < t < t k +i 
is obtained as well. Here W and B denotes two independent 1- 
dimensional Brownian motions. The increments Wt k+1 — Wt k and 
B tk+1 — B tk are simulated by quasi-random independent Gaussian vari- 
ables. 

The underlying price process of Heston model is given as follows; 

(dX t = rX t dt + ^V t X t dW t , 

\dV t = k(6 - V t )dt + uy/Vt(pdWt + v 7 ! - P 2 dB t ) 

for r, k, 9, v > and —1 < p < 1. Then the symmetrized path X is 
constructed clS EL solution to the following SDE; 



dt 



dX t = (rX t I {jtt>K} - r{2K - X t )I { x t<K} ) 

+ (Vv t X t I {jlt>K} - VVt(2K - X t )I {jtt<K} ) dW u 
where V is the solution of SDE (EDI). 
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The underlying asset price of A-SABR model is described as 

idX t = rX t dt + V t xfdW t , 

\dV t = X(0 - V t )dt + vV t {pdW t + y/1 - p 2 dB t ) 

for r, A, 0, v > 0, (3 > \ an d ~ 1 — P — 1- Then the symmetrized 
process X is given by the following SDE; 

d* t = (rX t J { ^ } - r(2K - X t )I {jtt<K} ) dt 

+ (v t xfl { x t>K} - V t (2K - X t fl {jtt<K} ) dW 1>t , 

where V is the solution of SDE (|2IJ1. 

Tables [7] - [10] below are simulation results of the price of a single 
barrier call option under Heston's and A-SABR model, respectively. 
We set the parameters as [X = 100, V = 0.03, K = 95, H = 
90, = 0.03, r = 0, T = 1, k = 1, p = -0.7, i/ = 0.03] in Heston 
model (Tableland Table ED, and [A" = 100, V = 0.5, S = 95, K = 
90, = 0.03, r = 0, T = 1, j3 = 0.75, A = 1.0, p = -0.7, = 0.3] in 
A-SABR model (Tableland Table ED, and 

Table [7]and[9j r = 0, 

Table HandE)) r = 0.02, 

in the experiments. Benchmark is given in the same setting of Section 
13.2.21 We again observe the superiority of our scheme. 



4.2 Application to Pricing Double Barrier Op- 
tions under the Stochastic Volatility Models 

Fix K, K' > 0. Let us consider a double barrier option knocked out 
if price process X hit either the boundary K or K + K' . The price 
of a double barrier option with payoff function / and barriers K and 
K + K' is given by E[f(X T )I {T(KK+K , )>T} ], where T( KiK+K i) is the 
first exit time of X from (K, K + K 1 ). In a similar way as the static 
hedging formula of a single barrier option, we obtain a static hedging 
formula if the price process satisfies APCS both at K and K + K'. 

Proposition 4.3 (jAkahori and Imamura 

(2Ql3)). IfX satisfies APCS 



at both K and K + K' , then for any bounded Borel function f and 
T > 0, we have 

v[f{x T )i {T{KK+Kl)>T} \ 

= ^ E[/(Xr - 2nK')I\ K+2 nK',K+(2n+l)K')( X T)] 

neZ ( 22 ) 

- ]T Bf(2K - (X T - 2nK'))I 

[K+(2n-l)K'K+2nK') 
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Table 7: Heston model 

X a = 100, V = 0.03, K = 95, H = 90, 9 = 0.03, r = 0, T = 1, k = 1, p = -0.7, k = 
0.03, benchmark of option price = 7.92706 



No. of 


No. of 






EM 


PCM 


simulation trials 


time steps 


EM 


PCM 


error (%) 


error (%) 


1000 


10 


8.638 


7.953 


9.0 


0.3 


8000 


20 


8.761 


8.167 


10.5 


3.0 


27000 


30 


8.466 


7.932 


6.8 


0.1 


64000 


40 


8.477 


8.017 


6.9 


1.1 


125000 


50 


8.366 


7.892 


5.5 


0.4 


216000 


60 


8.301 


7.877 


4.7 


0.6 


343000 


70 


8.246 


7.875 


4.0 


0.7 


512000 


80 


8.273 


7.902 


4.4 


0.3 


729000 


90 


8.221 


7.875 


3.7 


0.7 


1000000 


100 


8.212 


7.871 


3.6 


0.7 



Table 8: Heston model 

X = 100, V = 0.03, K = 95, H = 90, 6 = 0.03, r = 0.02, T = 1, k = 1, p = -0.7, v = 
0.03, benchmark of option price = 9.15602 



No. of 


No. of 






EM 


PCM 


simulation trials 


time steps 


EM 


PCM 


error (%) 


error (%) 


1000 


10 


10.308 


9.192 


12.6 


0.4 


8000 


20 


9.828 


9.197 


7.3 


0.5 


27000 


30 


9.572 


8.953 


4.5 


2.2 


64000 


40 


9.674 


9.133 


5.7 


0.3 


125000 


50 


9.632 


9.134 


5.2 


0.2 


216000 


60 


9.552 


9.093 


4.3 


0.7 


343000 


70 


9.525 


9.096 


4.0 


0.7 


512000 


80 


9.524 


9.135 


4.0 


0.2 


729000 


90 


9.498 


9.116 


3.7 


0.4 


1000000 


100 


9.454 


9.106 


3.2 


0.5 
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Table 9: A-SABR model 

X a = 100, V = 0.5, S = 95, K = 90, 6 = 0.03, r = 0, T = 1, [3 = 0.75, A = 1.0, 
—0.7, v = 0.3, benchmark of option price = 6.59534 



No. of 


No. of 






EM 


PCM 


simulation trials 


time steps 


EM 


PCM 


error (%) 


error (%) 


1000 


10 


6.643 


6.478 


0.7 


1.8 


8000 


20 


6.708 


6.591 


1.7 


0.1 


27000 


30 


6.701 


6.584 


1.6 


0.2 


64000 


40 


6.671 


6.565 


1.1 


0.5 


125000 


50 


6.668 


6.568 


1.1 


0.4 


216000 


60 


6.672 


6.581 


1.2 


0.2 


343000 


70 


6.669 


6.585 


1.1 


0.2 


512000 


80 


6.671 


6.597 


1.1 


0.0 


729000 


90 


6.655 


6.579 


0.9 


0.2 


1000000 


100 


6.646 


6.576 


0.8 


0.3 



Table 10: A-SABR model 

X = 100, V = 0.5, S = 95, K = 90, 9 = 0.03, r = 0.02, T = 1, /3 = 0.75, A = 1.0, 
—0.7, = 0.3, benchmark of option price = 8.71005 



No. of 
simulation trials 


No. of 
time steps 


EM 


PCM 


EM 
error (%) 


PCM 
error (%) 


1000 


10 


9.493 


8.779 


9.0 


0.8 


8000 


20 


9.081 


8.582 


4.3 


1.5 


27000 


30 


9.106 


8.723 


4.5 


0.2 


64000 


40 


9.029 


8.656 


3.7 


0.6 


125000 


50 


9.007 


8.683 


3.4 


0.3 


216000 


60 


9.008 


8.710 


3.4 


0.0 


343000 


70 


8.988 


8.707 


3.2 


0.0 


512000 


80 


8.940 


8.670 


2.6 


0.5 


729000 


90 


8.923 


8.671 


2.4 


0.4 


1000000 


100 


8.929 


8.680 


2.5 


0.3 
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Of the formula (|22p . the left-hand-side is the price of a barrier 
option, and the right-hand-side is an infinite series of the prices of 
plain-vanilla options. It means that a double barrier option can be 
hedged by infinite plain- vanilla options. Practically, the series should 
be approximated by finite terms. In our numerical scheme, however, 
finite sum approximation is not necessary as we will explain later in 
Remark 14.51 

We give a numerical scheme of a double barrier option under 
stochastic volatility model by using the symmetrized process which 
satisfies APCS at both K and K + K' . The scheme is summarized as 

Proposition 4.4. Set 

&n(x,v) 

= ^2 ail ( X ~ 2nK ' ' V ) I {K+2nK>,K+(2n+l)K>)(x) 

nez (23) 

- ^2 au ( 2K ~( x ~ 2nK ')i v ) I [K+(2n-l)K',K+2nK>){x), 

nez 

fa{x,v) 

= mix - 2nK', v)I [K+2nK > ,K+(2n+l)K')(x) 

nez (24) 

- XI - ( X - 2nK ')^)I[ K+( 2n-l)K',K+2nK'){x), 

nez 

and let X be the unique (weak) solution to 

dX t = & n (X t , V t )dW t + fii(X t , V t ) dt, 

where V is the solution to SDE \16\). Then, it holds for any bounded 
Borel function f and t > that 

E[f{X t )l {T{KK+Kl)>t] \ 

= yi E [f(Xt - 2 n K ')I\K+2nK>,K+(2n+l)K'){ X t)] 

- X E[f(2K - (X t - 2nK'))I 

[K+(2n-l)K' ,K+2nK') 

nez 

Proof This is an easy consequence of Proposition 14.31 □ 



Remark 4.5. The infinite series of the right hand side in I123\) and 
P^P is expressed by the following; 

(the right hand side of (|23[) ) 

f <T( x -[*jK.]K',v) ifi^^O mod2, 

I -a(2K-(x-([^]-l)K%v) if[^] = l mod 2, 
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and 

(the right hand side of 

li(x-[^]K',v) if[^} = mod 2, 

-^2K-(x-([^}-l)K'),v) if[^] = l mod 2. 

Therefore the discretized process of(X, V) by Euler-Maruyama scheme 
can be simulated without approximating the infinite series by finite 
sums. Similarly, we have 



(the right hand side of (|25|) ) 
= E f(X t \^l^]K\V t )I { ^ Q mod2} 

-f(2K - (X t - - 1)K% V t )I {{s ^ mod 2} 

Therefore Put-Call symmetry method is available for the pricing of a 
barrier option. 

Table [TT1 and Table [T2l below are numerical results of the pricing of 
a double barrier call option under Heston model and A-SABR model, 
respectively. We take 

Table HQ X = 100, V = 0.03, S = 95, K + K' = 115, K = 85, 9 = 
0.03, r = 0.02, T = 1, k = 1, p = -0.7, v = 0.03, 

Table H2JX = 100, V = 0.3, 5 = 95, if + K' = 110, if = 90, 9 = 
0.3, r = 0.02, T = 1, /3 = 0.75, A = 1, p = -0.7, i/ = 0.3, 

in the experiments. Benchmark is given by the same setting of Section 
13X21 

We still see that the put-call symmetry method beats the path-wise 
EM. 



5 Concluding Remark 

The new scheme, which is based on the symmetrization of diffusion 
process, is, though not theoretically, experimentally proven to be more 
effective than the path-wise Euler-Maruyama approximation scheme. 
The scheme is also applicable to stochastic volatility models including 
Heston's and SABR type. 
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Table 11: Heston model 

X Q = 100, V = 0.03, S = 95, K + K' = 115, K = 85, 6 = 0.03, r = 0.02, T = 1, k = 
1, p = —0.7, f = 0.03, benchmark of option price = 1.40319930 



No. of 
simulation trials 


No. of 
time steps 


EM 


PCM 


EM 
error (%) 


PCM 
error(%) 


1000 


10 


2 


987 


1.671 


112.869 


19.1 


8000 


20 


2 


368 


1.498 


68.759 


6.7 


27000 


30 


2 


144 


1.588 


52.785 


13.2 


64000 


40 


2 


045 


1.475 


45.770 


5.1 


125000 


50 


1 


921 


1.402 


36.903 


0.1 


216000 


60 


1 


876 


1.453 


33.662 


3.6 


343000 


70 


1 


820 


1.411 


29.728 


0.6 


512000 


80 


1 


792 


1.438 


27.733 


2.5 


729000 


90 


1 


765 


1.411 


25.791 


0.6 


1000000 


100 


1 


744 


1.416 


24.281 


0.9 



Table 12: A-SABR model 

X = 100, V = 0.3, S = 95, K + K' = 110, K = 90, 6 = 0.3, r = 0.02, T = 1, ft = 
0.75, A = 1, p = -0.7, v = 0.3, benchmark of option price = 2.46950606 



No. of 


No. of 






EM 


PCM 


simulation trials 


time steps 


EM 


PCM 


error (%) 


error (%) 


1000 


10 


3.779 


2.451 


53.017 


0.8 


8000 


20 


3.427 


2.566 


38.768 


3.9 


27000 


30 


3.164 


2.442 


28.126 


1.1 


64000 


40 


3.037 


2.489 


22.997 


0.8 


125000 


50 


2.955 


2.514 


19.640 


1.8 


216000 


60 


2.915 


2.480 


18.036 


0.4 


343000 


70 


2.875 


2.481 


16.438 


0.5 


512000 


80 


2.838 


2.478 


14.906 


0.4 


729000 


90 


2.806 


2.464 


13.631 


0.2 


1000000 


100 


2.779 


2.465 


12.540 


0.2 
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